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FOREWORD 

This Indian Standard was adopted by the Bureau of Indian Standards, after the draft finalized by the Rock 
Mechanics Sectional Committee had been approved by the Civil Engineering Division Council. 

Subsidence is a very comraan hazard in tunnelling works. This standard provides necessary information for 
predicting the subsidence in single seam coal mines. 

The composition of Committee responsible for the formulation of this standard is given at Annex A. 

For the purpose of deciding whether a particular requirement of this standard is complied with the final value, 
observed or calculated expressing the result of a test or analysis, shall be rounded off in accordance with 
IS 2 : 1960 'Rules for rounding off numerical values (revised)'. The number of significant places retained in the 
rounded off value should be the same as that of the specified value in this standard. 
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Indian Standard 

GUIDELINES FOR USE IN PREDICTION OF 

SUBSIDENCE AND ASSOCIATED PARAMETERS 

IN COAL MINES HAVING NEARLY HORIZONTAL 

SINGLE SEAM WORKINGS 



1 SCOPE 

1 . 1 This standard applies to the approximate prediction 
of subsidence in coal mines. Prediction of horizontal 
and vertical displacements, slope and horizontal 
strains for horizontal or nearly horizontal seams are 
dealt with in this standard. 

1.2 The provisions of this standard are applicable to a 
single seam mine at a time. 



2 TERMINOLOGY 

2.0 For the purpose of this standard the following 
definitions shall apply. 

2.1 Apparent Angle of Dip (a) — The apparent angle 
of dip of the seam with respect to horizontal plane. 

2.2 Angle of Draw (^ — The angle between the 
vertical line and the line connecting the upper edge of 
the working to the edge of the subsidence trough 
(^eeFig. 1). 
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2.3 Subsidence (5) — The vertical movement of any 
point on the surface caused by an underground 
excavation (see Fig. 1). 

2.4 Maximum Subsidence {SJ — The maximum 
vertical displacement caused by an underground 
excavation on the surface in a subsidence affected area 

(see Fig. 1). 

2.5 Maximum Possible Subsidence (S^^^) — The 
maximum vertical movement that can be caused on 
the surface in a given locality by an underground 
excavation of such a critical dimension (see Fig. 1). 

2.6 Panel — A part of the mine workings from where 
the mineral or coal has been or is being mined out. 

2.7 Pillar — A block of coal left unexcavated to offer 
support to the roof of the workings. 

2.8 Critical Length or Width — The length or width 
of excavation or working panel that causes maximum 
possible subsidence (5^;^^) at a single point on the 
surface (see Fig. 1). 

2.9 Super-Critical Length or Width — The length 
or width of the panel larger than critical dimension 
(see Fig. 1). 

2.10 Sub-Critical Length or Width — The length 
or width of the panel smaller than critical dimension 
(see Fig. 1). 

2.1 1 Curvature (l/p) — The reciprocal of the radius 
of curvature of any part of the subsided surface profile. 

2.12 Depth (d) of Overburden — The vertical 
distance from the roof of the extracted seam to a point 
vertically above it on the surface. 

2.13 Slope (g) — The slope of any part of the 
subsidence profile or subsidence trough. 

2.14 Strain (ej — Change per unit length in the 
distance between any two points on the surface in a 
specified direction. The tensile strains are considered 
as positive (ej and the compressive strains as negative 
(ej. 

2.15 Critical Area (A^^.J — The minimum area of 
excavation that causes maximum possible subsidence 
of one point on the surface (see Fig. 1). 

2.16 Sub-Critical Area (A^^^) — An excavation area 
smaller than the critical area. The subsidence caused 
by this area is smaller than maximum possible 
subsidence (Sj^^^) at any point (see Fig. 1). 

2.17 Super-Critical Area (^.^p) — An area of 
excavation larger than the critical area. The subsidence 
caused by this area of excavation is equal to the 



maximum possible subsidence (iS^^^,^) at floor of the 
subsidence trough (see Fig. 1). 

3 PREDICTION OF MAXIMUM POSSIBLE 
SUBSIDENCE (S^J 

3.1 Data Required for the Prediction of 5^^^ 

The following data is required for calculating the 
anticipated magnitude of 5^^^ for a given mining and 
geological situation: 

a) Extracted thickness of the seam (m); 

b) Extraction ratio (ER) that is the ratio of 
extracted volume of coal or mineral to its total 
volume in the given panel; 

c) Method of goaf support that is caving, 
stowing, etc; 

d) Depth of workings (d) below the surface leve I ; 

e) Angle of dip (a) of the seam; and 

f) The Young's modulus £^(in the direction 
perpendicular to the bedding plane) and the 
shear modulus G, of the rockmass. In the 
absence of homogeneity; the overall 
(weighted average) values ofE^ and G should 
be used for the following situations: 

1) Presence of hard rock layers, namely, 
sandstone, limestone, etc, in the 
overburden rockmass; 

2) Presence of old or current workings in 
the vicinity of the panel and the parting 
between them; and 

3) Presence of geological discontinuities, 
namely, joints, fractures, etc. 

3.2 Prediction of 5'J^^^^ 

The maximum possible subsidence in a given area 
due to critical or super-critical excavation is given 
by: 



^Ma.^ f"^fSfRAd't 



(1) 



where 
m 



- Extracted seam thickness, 

= Extraction factor, 

= Goaf treatment factor, 

= Rock factor for the combined effect of 

composition and condition of the 

overburden rockmass, 
= Factor for the effect of depth of working, 
= Factor for the effect of dip of the seam, 

and 
= Time factor (taken to be unity for finished 

subsidence). 
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3.2.1 Extraction Factor {e^ 

The effect of partial extraction on S^^^ is given by the 
extraction factor equal to {ER)^ where ER is the 
extraction ratio (ratio of extracted volume to the total 
volume of the coal within the panel, in fraction) and 
A" is a constant. The value of /^ may be taken equal to 
1 for general cases of soft coal seam and 2 for hard 
coal seams and may be interpolated between 1 and 2 
for intermediate cases. Alternatively, the value of 
extraction factor should be taken from Fig. 2. Coal 
having uniaxial compressive strength of 15 MPa 
( 1 50 kg/cm^) or more should be considered as of 'hard' 
quality and less than that of 'soft'. 

3.2.2 Goaf Treatment Factor (g^) 

The value of this factor should be taken as 0.95 for 
caving cases and 0.07 to 0.10 for cases of hydraulic 
sand-stowing as goaf support. 

3.2.3 Depth Factor (d^) 

The factor for the depth of coal seam should be taken 
equal to 0.87 for depths up to 250 m; 0,96 for depths 
between 250 and 400 m and 1 ,0 for depths beyond 400 m. 

3.2.4 Factor for the Apparent Dip of the Seam {d ) 

This factor should be taken as equal to cosa, where a 
is the angle of apparent dip in degrees (for its values 
up to 20"). 

3.2.5 Rock Factor (R^) 

The effect of composition and condition of overburden 
rockmass on S^^^ should be quantified by a factor called 



rock factor R^. whose value lies between for no 
subsidence to 1 ,0 for the maximum value ofS,^ . The 

Max 

presence of hard rock layers in the overburden 
rockmass results in the reduction of subsidence, 
whereas its fragmentation makes it weaker, thus 
increasing the subsidence. 

On the basis of its condition, the overburden rockmass 
is classified from subsidence point of view into five 
categories as given in Table 1. In Table 1, £^ is the 
Young's modulus of overburden rockmass in a 
direction perpendicular to the bedding plane and G is 
shear modulus. The ratio G/E^ represents the strength 
of rockmass in which E^ may be empirically corelated 
to the presence of hard rock layers in the overburden 
by the following formula: 

E^(\n MPa) = 500 + (195 x percentage of hard 
rock layers in the over 
burden) 

The condition of the overburden rockmass that is the 
presence of natural discontinuities and the degree of 
fragmentation caused in it because of repeated working 
may be expressed in terms of G/E^ ratio. The range of 
G/E^ ratio for each classification of rockmass is also 
given. The rock classification for a given overburden 
rockmass should be found out either from its G/E 

V 

ratio or from the description given in Table 1 . 

The rock factor R^ should be obtained for a given 
rockmass (as classified in Table 1), either by knowing 
its Young's modulus £^and shear modulus G, or by 
knowing the percentage of hard rock layers in the 
overburden fi^om the respective curve given in Fig. 3, 
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Fig. 2 Effect of Partial Extraction on Subsidence 
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Fig. 3 Relationship Between R^ and Hard Rock Layers in the Overburden 



Table 1 Classification of Overburden Rockmass 

{Clause 32.5) 



SI 


Classification 


G/E, 


Description 


No. 








(1) 


(2) 


(3) 


(4) 


i) 


Competent 


0.07-0.38 


Massive rockmass with few 




rockmass 




planes of weaknesses 


ii) 


Undisturbed 


0.03-0.07 


No previous mining but few 




rockmass 




natural discontinuities 


iii) 


Partially 


0.016-0.03 


Parting between the mined 




disturbed 




seams more than 5 times the 




rockmass 




seam thickness; or cases 
with no previous mining but 
many natural discontinuities 


iv) 


Disturbed 


0.005- 


Parting between the mined 




rockmass 


0.016 


seams not more than 5 
times the seam thickness; or 
having over lying/under 
lying or side lying workings 
in the vicinity of the panel. 


V) 


Highly disturbed 


< 0.005 


Highly fragmented 




rockmass 




rockmass witli repetitive 
workings or having very 
thick seams mined in 
descending slicing 



4 PREDICTION OF MAXIMUM SUBSIDENCE 
FOR A SUB-CRITICAL AREA OF EXTRACTION 

4.0 The maximum possible subsidence S^^ occurs at 
or near the centre of the panel when its length (0 and 
width (PF), each is more than or equal to critical 
dimension. 

4.1 If the length of the panel is critical or super-critical 
and width is sub-critical, the maximum subsidence S^ 
near the centre of the panel will be less than S^^^ and 
is given by: 

(2) 

(3) 



where 



s-s^w 



W' = [l- exp{-n(f^^df}] 



4.2 If besides the width, the length of the 
panel is also of sub-critical dimension, the maximum 
subsidence 5„ will be given by: 



where 



S =S,, W u 

Max 



L' = [l- exp{ -n(l/df}] 



(4) 
(5) 
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where n is an empirical constant value of which lies 
between 2.5 and 3.5 (average 3.0) 

5 SUBSIDENCE PROFILE FOR A 
RECTANGULAR PANEL 

The subsidence at any point along the central line in 
a rectangular panel and along the other lines parallel 
or perpendicular to the central line can be predicted 
by knowing: 

a) the maximum subsidence due to the 
excavated panel, 

b) the maximum subsidence along the given 
line, 

c) distance of the line from the point undergoing 
maximum subsidence in whole of the panel, 

d) the distance of the given point from the point 
undergoing maximum subsidence along that 
line, and 

e) the angle of draw on both sides of the 
subsidence profile. 

5.1 Angle of Draw (^ 

Figure 4 gives the magnitude of angle of draw ^, in 
degrees on the commencing side (or static end) for 
different W/d ratios and for a given rockmass 
represented by rock factor R^ The angle of draw ^^ on 



the other end (dynamic end or moving face side) 
should be found out from Fig. 5. 

5.2 Subsidence Trough Due to Nearly Horizontal 
Rectangular Excavation 

The subsidence trough in a horizontal or nearly 
horizontal excavation goes upto a distance r^ equal to 
cf tan 1^, on the commencing side and up to a distance 
r^ equal to c/ tan ^^ beyond the limit of excavation on 
the dynamic face side as shown in Fig. 6. 

The point undergoing maximum subsidence S^ or 5"^,^^^ 
lies in the centre of a sub-critical and critical panel 
respectively. In case of super-ciritical panel, an area 
on the surface above the central portion of the panel 
undergoes the maximum possible susbsidence. 
Therefore in such cases a subsidence trough with a 
level floor at full subsidence S,^^^ is produced. The 
points undergoing full subsidence along a longitudinal 
or transverse central line profile are situated between 
two points^ ' and A"' {see Fig. 6) lying at a distance 
/•, equal to d tan ^, and at a distance r^ equal to 
dim ^2 respectively from the stationary rib and moving 
towards the centre of the panel. Alternately these two 
points X 'and X ' are situated at a distance of ( WI2 ~r^) 
and {WI2 - r^ from centre of the panel towards the 
commencing and finishing sides respectively {see Fig. 
6), where w is the width of the panel. 



Rf (0-91- 100) o 



UJ 

o 



o 



^0 



-z. 

UJ 

S. 30 
O u) 

5 Q 20 

OMT 

o 

UJ 10 - 

_i 

o 

z 
< 




1 



2 3 

W/d RATIO 



Fig. 4 Angle of Draw on Commencing Side (^) for Different WId Ratios and R^ Values 

5 



IS 15180 : 2002 



13 
Z 

■3Z 

to 



AOr 



30 






U. UJ 

o9 

_J 

o 
z 
< 



10 - 



Rf (0-91 -1-00) 




Rf (0 66 -0-9) 



o o 



Rf (0-36-0 -65) 



Rf (0-0-35) 



■35' 



28" 



-flL20 



.10" 



1 2 3 A 5 6 

W/d RATIO 

Fig. 5 Angle gf Draw on Finishing Side (^ ) for Different fV/d Ratios and /?, Values 



2 



W 



[_r, -^(^r,^ 




Fig. 6 Subsidence Trough Due to Super- critical Excavation 



5.3 Subsidence Profile for Rectangular Excavation 
in Nearly Horizontal Seams 

The vertical displacements 5"^ that is subsidence at a 
point at a distance x from the centre of a sub-critical 
panel may be found as shown in Fig. 8: 

5 =5 [exp{-M(x/(r + x)M] -.(6) 

and, for critical and super-critical extractions; 



where 



S=S^lsxp { -Mix/ir + xy}] . . .(7) 



distance of the given point from the centre 
of a critical or sub-critical panel, or 
distance of the given point from the nearest 
point undergoing S^^^ along that profile 



line (^see Fig. 6) in case of super-critical 
panel. 

r = Critical diameter = (/•,+ r^ 

M = Profileconstant on which the shape of the 
subsidence profile depends. Value of M 
may be found out from Fig. 7 depending 
on the value ofR^. The higher value of R^ 
(that is, for more disturbed rockmass) the 
more is the value of A/and the sides of the 
subsidence trough will be steeper. 

5.4 Subsidence Trough for Panels with Irregular 
Boundaries 

An example of extraction area with irregular boundary 
is shown in Fig. 8. The value of 5^^^ or S^ should be 
found out from equations (1) and (4) by averaging the 
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O Si SUBSIDENCE AT THE CENTRES OF IMAGINARY PROFILES 
• Sj SUBSIDENCE AT THE CENTRES OF ACTUAL PROFILES 

imaginary profiles considering rectangular boundary 

actual profiles considering actual boundary 

Fig. 8 Imaginary and Actual Profile Lines in Excavations with Irregular Boundary 



boundaries of the excavation and simulating it to be 
of rectangular shape. Two main profile lines one 
longitudinal and the other at right angles to this and 
both passing through the centre of the panel should 
be drawn. 

The extraction area should be treated as rectangular 
in shape, considering the longest dimension as its 
length and the longest dimension perpendicular to this 
as its width. The value of subsidence S^ for all points 
at a suitable distance interval dx along these two main 
profiles are calculated by profile function that is from 
equations (6) and (7). Imaginary profile lines, parallel 
and equal in length to the main transverse profile are 
drawn (at a distance dx apart) and undergoing 



subsidence S^ at their centres are calculated. The values 
of subsidence S. at all the points at a distance dx apart 
along these imaginary profiles (with respect to the 
assumed rectangular extraction) are calculated with 
equations (6) and (7) given in 5.3 by replacing S.^ for 
S or S.. . The subsidence values S. thus obtained are 

Max J 

due to assumed rectangular excavation. 

Now considering each transverse profile of its actual 
length (or width) intercepted between actual 
excavation boundary, their centres are found out at 
which subsidence value S- is known already. The 
magnitude of angle of draw for each actual width (or 
length) ratio is found out from Fig. 4 and Fig. 5 for 
calculating the extent of each profile. The subsidence 
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Fig. 9 Subsidence Profiles by Interpolation 



at all points at a distance dx apart along these actual 
profiles {see Fig. 9) is now estimated from equations 

(6) and (7) by substituting S. for 5^^^. 

6 PREDICTION OF OTHER ASSOCIATED 
PARAMETERS 

6.1 The parameters associated with subsidence 
movements, which may have damaging impacts on 
various surface properties and structures are, the slope 
of the subsidence profile curve, horizontal 
displacements of the surface points, horizontal strains 
the ground points experience and the curvature of the 
ground. 

6.2 Determination of Slope or Tilt 

Slope or tilt gives the variation of differential 
settlement between two points on the ground. The 
slope g., of a subsidence profile at the rth point at a 
distance x from the centre of the panel is the derivative 
of subsidence with respect to the distance, ITie slope 
between two consecutive points say /th and (/ - l)th 
from the panel centre can geometrically be obtained 
as; 

g.^M ^ ^ {s-s, ,)/A (8) 

where 

s. and i;_, are the vertical settlements (or subsidence) 
of the /th and (/' - l)th points and A^^ is the distance 
between these two points. Therefore the slope is the 
differential subsidence between two consecutive points 
on subsidence profile situated at 1 m apart and is 
expressed in mm/m. 

6.3 Determination of Horizontal Displacement 

The surface points in a subsidence trough are subjected 



to vertical as well as horizontal displacements. The 
profile curve of horizontal displacements is similar in 
nature to the profile curve of slopes. A linear 
proportionality exists between the two curves. The 
horizontal displacement u. at rth point from the panel 
centre is given by 

«i = figi (9) 



where 
B 



slope at the rth point, 

Proportionality constant which is 
dependent on the nature of the overburden 
rockmass that is the rock factor R^ and on 
the W/d ratio of the extracted panel {see 
Fig. 10). 



6.4 Determination of Horizontal Strains 

The horizontal strains {e.j are given by: 

Am, 



^i±- 



(10) 



where 



Difference in horizontal displacements of 
the two consecutive points that is rth and 
(/- l)th which were originally situated at 
a distance A^ apart. Am. is found out as 
follows: 



Am. = u-u, 



i-l 



.(11) 



The horizontal strains as obtained from equation ( 1 0) 
are positive (ej or tensile for points situated in the 
outer areas of the subsidence trough and are negative 
or compressive (ej for points situated towards the 
panel centre. 
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Fig. 1 Values of B for Different W/d Ratio 



6.5 Determination of Surface Curvature 

The curvature of the ground {K^ is the first derivative 
of the slope. It is the change in the slope per unit 
length. 

The curvature is given by. 



K. 



' A. 



.(12) 



where 

g. = difference in slope at rth and (/ - l)th 
point. 

7 NON-EFFECTIVE WIDTH 

The maximum width of an underground opening that 



does not produce any significant movement on the 
ground surface is called the non-effective width 
(NEW). TTie concept of NEW is helpful in the design 
of the underground excavations and in controlling or 
minimizing the subsidence movements by keeping 
them within tolerable limits. The value of NEW is 
normally expressed as the ratio of this width to the 
depth of the excavation. The value of NEW depends 
on the nature of overburden rockmass, the thickness 
of the coal seam and the extraction ratio. The values 
of NEW for longwall workings (100 percent 
extraction) and for 50 percent extraction in board and 
pillar workmgs are given in Fig. 11 and Fig. 12 for 
indifferent seam thicknesses. 
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NEW — NON EFFECTIVE WIDTH 
NEW FOR 50 % EXTRACTION IN BOARD 
AND PILLAR WORKINGS 
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Fig. 1 1 Values of NEW for 50 Percent Extraction 
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